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Abstract 



We investigate the magnetic susceptibility of the QCD vacuum, based on the instanton vacuum. 
Starting from the instanton liquid model for the instanton vacuum, we derive the light-quark 
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partition function Z[V,T,rh] in the presence of the current quark mass rh as well as the external 
Abelian vector and tensor fields. We calculate a two-point correlation function relevant for the 
magnetic susceptibility and derive it beyond the chiral limit. We obtain for the different flavors the 
following magnetic susceptibility: Xu,d{i^ u d^Pu^o ~ 40 ~ 45MeV, while X s (m/Ws)o — 6 ~ lOMeV 
j> ' with the quark condensate (ix/}'iJj)o. 
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1. The QCD vacuum is one of the most complicated states filled not only by perturbative 
but also strong nonperturbative fluctuations. In particular, the condensates of the quarks 
and gluons characterize its nonperturbative aspect. The quark condensate plays a role of 
an order parameter associated with spontaneous chiral symmetry breaking (S^SB) which is 
probably one of the most important quantities for low-energy hadronic phenomena. 

The presence of a constant external electromagnetic field induces another type of 
condensates leading to the nonzero magnetic susceptibility Xf °f the QCD vacuum defined 
as: 

(0|V>}cvV/|0)f = e /Xf WfTl>f)o F i»» (!) 

where e/ denotes the quark electric charge, / is the corresponding flavor. It is natural 
to have the quark condensate (iipfipf)o in the chiral limit as a normalization factor in the 
right-hand side of Eq. (JTJ), since S^SB is responsible for these quantities, i.e. the quark 
condensate and magnetic susceptibility. In a recent paper [jj it was suggested that the 
magnetic susceptibility Xf ma Y be measured in the exclusive photoproduction of hard dijets 
7 + iV — > (qq) + N. The value of Xf was predicted by using the method of the QCD sum 
rule and vector dominance 

BIB: XfWtyfh = 40 - 70 MeV at the scale of 1 GeV. 

In the present work, we want to investigate the magnetic susceptibility Xf °f the QCD 
vacuum within the framework of the instanton vacuum. Since the instanton vacuum explains 
S^SB naturally via quark zero modes, it may provide a good framework to study the Xf- 
Moreover, there are only two parameters in this approach, namely the average instanton size 
p « g fm and average inter-instanton distance R ~ 1 fm. The normalization scale of this 
approach can be defined by the average size of instantons and is approximately equal to p _1 w 
0.6 GeV. The values of the p and R were estimated many years ago phenomenologically by 
Ref. as well as theoretically by Ref. 0. Furthermore, it was recently confirmed in various 
lattice simulations of the QCD vacuum p], 0, 0] . Very recent lattice calculations of the quark 
propagator 0, 11 1 are in a remarkable agreement with that of Ref. 

The magnetic susceptibility of the QCD vacuum was estimated already in Ref. |2j 
however, done in the chiral limit and without the consideration of the current conservation. 
Since we want to investigate it beyond the chiral limit, we first need to extend Refs. 0, 
Eii . I3 | and Lee-Bardeen formula ^| for a low-frequency (p <C p^ 1 ) part of the quark 



determinant: 

Deti ow = detS, B {j = imS^ + a jh a_ + = -($_ |^|$+,o), (2) 

taking into account the current quark mass. Here, m denotes a current quark mass, which is 
assumed to be small, a^- stands for the overlapping matrix element of the quark zero modes 
$±,0 generated by instantons (antiinstantons). In general, we define $±, n by the following 
equation: 

(ip + fi±)Q±, n = \ n ®±, n , (3) 

where ^4 M) ± is an instanton (antiinstanton) field and Ao = by convention. The matrix 
element is nonzero only between instantons and antiinstantons and vice versa due to 
specific chiral properties of the zero modes. 

The overlapping of the quark zero modes allows quarks to propagate through the 
instanton-antiinstanton medium by jumping from one instanton to another one. Hence, 
Deti ow is reduced to the determinant of the finite matrix B in the subspace of zero modes. 
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This result was reproduced by Refs. jl3l . 16] in different ways. It is well known that the 
fermionic determinant is proportional to m) N+ ~ N -\ and strongly suppresses the fluctuations 
of \N + — iV_|. Therefore, we can take N + = AL = N/2 in the final formulas. 

In order to calculate the magnetic susceptibility, we have to find the light-quark partition 
function Z[V,T,rh] in the presence of the Abelian external vector and antisymmetric 
tensor fields. So, we first calculate a zero-mode approximated quark propagator S in 
the instanton ensemble and in the presence of these external fields. Using this propagator, 
we are able to compute a low-frequency part of the quark determinant Det[ ow and arrive at 
the extension of Eq. (j2J) for the non-zero m with the external Abelian vector V and tensor 
T fields. 

The smallness of the packing parameter vr(-^) 4 « 0.1 makes it possible to average the 
determinant over collective coordinates of instantons with fermionic quasiparticles, i.e. con- 
stituent quarks ip introduced. The averaged determinant turns out to be the light-quark 
partition function Z[V, T, m] which is a functional of V and T and can be represented by a 
functional integral over the constituent quark fields with the gauged effective chiral action 
S[if>\ if), V, T]. However, it is not trivial to make the action gauge-invariant due to the non- 
locality of the quark-quark interactions generated by instantons. In the previous paper 0| , 
it was demonstrated how to gauge the nonlocal effective chiral action in the presence of 
the external electromagnetic field. The gauged action was successfully tested by showing 
that the low-energy theorem of the axial anomaly related to the process GG — > 77 (see also 
0G3) is satisfied. 

In the present work we refine the derivation of the gauged effective chiral action and 
apply it to the calculations of the magnetic susceptibility Xf °f the QCD vacuum, taking 
into account both current quark mass and nonlocal currents. 

2. It is clear that the magnetic susceptibility Xf is related to the two-point correlation 
function of the quark tensor and vector currents at the zero momentum. In order to calculate 
this correlation function, we need to treat the light-quark partition function Z[V,T,rh] as 
a functional of external abelian vector and tensor T^ v fields. The relevant total quark 
determinant 

Det := Det(z0 + A + ef + a^T^ + irh) (4) 

can be splitted into two parts which correspond to low and high momenta, respectivley, 
with respect to some auxiliary parameter M\ lying inside an interval, -C M\ <C p^ 1 : 
Det = Detiow x Dethigh [12]. Dethigh comes from fermion modes with Dirac eigenvalues from 
the interval Mi to the Pauli-Villars mass M, and Deti ow accounts for eigenvalues less than 
Mi. The product of these determinants is independent of the scale Mi. However, we can 
deal with them only approximately. The high-momentum part Dethigh can be written as a 
product of the determinants in the field of individual instantons, while the low-momentum 
one Det^w has to be treated approximately, would-be_ zero modes being only taken into 
account. In Ref. fl2l |. it was demonstrated that the Deti ow depends weakly on Mi in the 
wide range, which serves as a check for this approximation. Based on this prescription, we 
now calculate the low-momentum part of the quark determinant in the presence of instantons 
as well as external V^, T^ u fields, and then average it over the collective coordinates of the 
instantons in order to find the effective low-energy QCD partition function as a functional 
of the external field Z [V, T, rh] . 

We first define the total quark propagator S in the presence of the instanton ensemble A 
and external fields V, T and the quark propagator Si with a single instanton Aj as well as 
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V and T: 

S ~ ift + gA + ef + a ■ T + im Si ~ ift + gfii + ef + a ■ T + im ^ 

where a ■ T = o^JT^. We assume that the total instanton field A may be approximated 
as a sum of the single instanton fields, A = Y2i=i ^> which is justified with the above- 
mentioned average size of instantons p ~ 1/3 fm and average inter-instanton distance R ~ 
1 fm. Defining the quark propagator So with external fields V and T and the free one So as 
follows: 

1 1 

^° = — v7~ 7ri~, ~i S = rx-, (6) 

ip + ey + a ■ 1 + im ip + %m 

we can expand the quark propagator S with respect to a single instanton: 

S = S + J^iSi - So) + - S )S \S 3 - So) + ■ ■ ■ ■ (7) 

In order to specify the gauge dependence, we rewrite Si and So in the following form: 

Si = LiS i L i , Sj = ~rz ■ t: ■ rrr — =— —j 

2jZ7 + p/lj + ey' + a - T + im 

So = L i S' L i x , = t^— — — — , (8) 

ip + ey + a ■ 1 + im 

where ef' = L~ 1 (i^ + ey r )L and the gauge connection L can be written as the path-ordered 
exponent: 

Li(x, zt) = P exp (ie jf ^^(Oj , (9) 

where Zj denotes an instanton position. 

In the case of the small quark mass m, we can expand the quark propagator Si 2| as 
follows: 

Si(m >• 0) = tit H ~ . (10) 

ip im 

While it gives a proper value for the (&oi\Si(m — > 0)|$oi) = the second term has wrong 
chiral properties: 

5i(^_ > 0)|$ w > = ^ + i|««>. (11) 

However, we may neglect this term as m — > 0. 

If the m is not small, then we propose the following approximation: 

S t = So + Soip ^^ ipSo, (12) 
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where 

d = -($oi\ipS ip\$oi) = im(^oi\S ip\^oi}. (13) 
The merit of the approximation proposed in Eq. (|12|) lies in projecting Si onto the zero-mode: 

= i|$oi), ($ i|^ = ($ 0i |tV. (14) 

im im 

Including the external fields V and T, we get 

Si = S 0i + SoiLiip -^—^- ipL^Soi, (15) 

Q Oi 

where 

{^i\ip{Lr l S Q Li)ip\%^. (16) 
Thus, we can express Eq.flTJ) as 

S = So + SaJ^Liimio) (z^Tr) (^l^^o, (17) 

where 

Aj + 2ij = ($ 0i \ip(Lr 1 S L j )ip\$ oj ), D tj = {h - Ci )<%. (18) 
Introducing the following function: 

which has the same chiral properties as the zero-mode eigenfunction |$oi) ; and taking the 
trace over the subspace of the zero-modes, we obtain the following expression: 

fr(5-5 o ) = -^(0o, J |^(4 2 )^l^)(0o,d(TTi^)l0o, J ), (20) 
~J ipS ip 

where Tr denotes the trace over the zero-mode subspace only. Bringing now in the matrix 

B(m)ij = (<poMipSoip)\M, (21) 
we get 

~ f m ~ ~ /"^ (m) 1 B(m) 

Tr / idm'(S(m') - S {m')) = Tr / dB{m')-„ = Trln . 1 ; , (22) 

J Mi Jb{Mi) B{rri) B(Mi) 

from which we find that B = B(m) is equivalent to a mere extension of the Lee-Bardeen 
matrix B in the presence of the external Abelian fields V, T and current quark mass 
rh. Thus, we derive for arbitrary flavor Nf 

M = (%S@L\ l S of L fj zp\% J ) (23) 



S °f = WT Tfr-—, L fi = Pexp ie f / d£MO . (24) 

ip + e f y + a - T + im f \ J Zi J 

If we introduce Grassmanian variables Vt{ and f^, then we can express the fermionic deter- 
minant of B as a functional integral over fermion fields: 

det£ = J dttdttexp (OBfi), (OBft)] = fif {^ipS^jip^)^ 

= (Det(ip + ef + a-T + im)y 1 jY[ d£l{ dVt{ D?P f Dip ] f 

ij 

J dx{rf(x) {4+e f f+ar-T + im f )rl>f (x)+fj{ / (x) + ^ (as) r,( (as) } ^ ^ j 



X e 



where the source fields 77/ and 77J are defined as: 



'// ^/(oo, ip, ri = ip\M^. (26) 

The integration over the Grassmanian variables Q and fl provides finally the low-frequency 
part of the quark determinant in the following form: 

Detiow = det B = (det(ip + ef + a ■ T + im))' 1 I JJ Dip f Dip j (27) 

J f 

( N+ N- 



x e 



where 

V ±J [^ f ^f} = J d 4 x {f)(x) Lf(x,z)ip$±,v(xrt ± )) j d 4 y (4 )0 ( 2 /;e ± )(^L+( 2 /^)V/(2/))28) 

It is obvious that V±j[Tp\ip] indicates the nonlocal interaction between constituent quarks 
generated by instantons. Since the range of the integration in Eq.(J2BJ is truncated at p, 
which is defined by zero- mode functions $±,0? the range of the nonlocality is determined by 
P- 

Note that the external vector field gauges not only the kinetic term of the effective 
action but also its interaction one reason is clear: Since the 

instanton-induced interaction is nonlocal, the gauge connection Lf must be attached to 
each fermionic line so that the effective action may be gauge-invariant. Unfortunately, we 
encounter an ambiguity arising from this gauge connection, since it is path-dependent in 
general. However, we will show now that the most preferable path can be found in the 
present case. 

3. In order to solve the problem of gauge dependence due to the nonlocal interaction 
induced by instantons, we have employed the gauge connection Lf. Though it causes in 
general some arbitrariness due to its path dependence, we can prove that such dependence 
can be minimized in the present work. To be more specific, we consider the extended zero 
mode: 

(ip + f)|<£ > = 0, ip = ip + /L, \^) = \%)-S NZ y\%), (29) 
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where A denotes an instanton field located at z, ) stands for the solution to order 
0(V), and Snz the well-known non-zero mode of the propagator in the instanton field 
(see the review [5] and references therein). Putting the gauge connection L(x, z) defined in 
Eq. © into Eq. J22), we obtain: 

(ip + ni*o> = °. V^L-^fy + VjL, \tf ) = L- 1 (x,z)\$ ). (30) 

ip |$ ) = and S NZ ip = 1 - |<E> )($o| being used, the solution |$q ) to order 0(V) in 
Eq. (jBlIj) can be reduced to the corresponding solution \&^) to order 0{V) without any 
trace of the path dependence arising from the gauge connection L. 

However, if the zero-mode approximation S^z ~ ■h — S 00 is used to find the solution to 

id 

order 0(V) 

I^SVW-W, (31) 

then applying the inverse gauge connection L^ 1 to Eq. (j3*T|) leads to 

|$$) = |$ ) - 5 of|$o) + Soo^ J d£M0)\$ ), (32) 

which depends on the path we choose explicitly. Thus, it is the best way that we choose the 
path minimizing — $go' ) || 2 . Expressing the external vector field in a Fourier form: 



V,(x) = j d*qV„(q)e^, (33) 

we can parameterize the path from x to z as follows: 

£p(x, z, t) = Zf, + (x - z)^t + p^x, z, t) (34) 

with parameter t G (0, 1) and boundary conditions p^(t = 0) = p^(t = 1) = 0. It is the 
most optimized path we can take, since p^ should depend only on the path difference x — z 
in order to satisfy the conservation of the linear momenta. More explicitly, we write down 
the path in the following form: 

£n(x, z,t) = z^ + (x - z) tl t+p fl t(l - t). (35) 

Then we have 

i j <%M0 = i j (PqV^e** j dt[{x - z), + (1 - 2t)^] e i <H(*-*)*+^(l-t)]. (36 ) 

In principle, we can express Eq. f!36|) in terms of the error function. However, we shall 
consider a simple case of qp <C 1 so that we can expand the vector current with respect to 
(q ■ p) and q • (x — z): 



Vpix) = / d^qV^qY^ ~ / ^qV^qy^ 



1 + iq ■ (x — z) — 



(q-{x- z)) 



21 



(37) 
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Substituting this expansion into Eqs. (j2Sj) and (|32*j) and carrying out some tedious but 
straightforward calculations, we obtain 

0o ] -^\\ 2 *Cp*N f J d^V\ q ){q-pf + --- (38) 

where the constant C = +0.0063 > and higher-order terms are neglected. Therefore, we 
can conclude that the minimum is achieved by taking the straight-line path that will be 
employed in the sequel. 

Moreover, since we are only interested in the case of a constant electromagnetic field 
strength F pu , we can even show that the magnetic susceptibility does not depend at all on 
the path. Let us consider the following difference: 

S (jf W(o) := J c ^M0 ~ [ dZMO, (39) 

where f c d^V^) is the line integral along the path C parameterized by £ M (x, z, t) = z p + 
(x—z) t j,t+a f j,(x—z, t) and J s _^ d^^V^) is calculated along the straight-line. This difference is 
reducible to the integral over the surface between these two contours due to Stoke's theorem: 



5 Uz^ M ®) := f d ^ V ^) = j dS^F^) = F^J d^ 



where the integral J c^a^) depends only on a distance (x— z) but not on x and z separately 
due to translational invariance. Denoting G pv {x — z) = f d^^a v {^)F{x — z), where F(x — z) 
represents the form factor F in coordinate space, we find that the change of the magnetic 
susceptibility 5xf due to the deformation of the contour turns out to vanish: 

F Xp Ti[S(p)G pX (p)F(p)S(p)a^ + S(p)F(p)G pX (p)S(p)a llu } = 0. (41) 

Thus, we are now allowed to take any path without changing the final results of the magnetic 
susceptibility. 

4. We are now in a position to derive the relevant partition function for the magnetic sus- 
ceptibility of the QCD vacuum. We first have to average the low-frequency part of the quark 
determinant deti ow over collective coordinates £±. Having exponentiated and bosonized the 
interaction terms [lil . Il3l 17 . liR 12(3 ]. we arrive at the effective chiral action with external 



fields pertinent to the magnetic susceptibility. Since the density of the instanton medium 
(it 2 (-|) fa 0.1) is low, we are able to average over positions and orientations of the instan- 
tons independently. In the next step we have to introduce auxiliary fields, i.e. the couplings 
A± for exponentiation and meson fields <&±(x) for bosonization, respectively. Having in- 
tegrated over fermion fields, we end up with the partition function in the presence of the 
external vector and tensor fields: 

Z[V,T,m} = J d\ + d\^D$ + D$^exp(-Sl\±,®±}) , (42) 
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where 
S[X±, 



iVj. In 



iV f 



iV 4 



V\ 



— N. 



± 



± 



(43) 



f d 4 xJ2( N f - l)X± Nf ' 1 (det^ ± (x)) N f- 1 



-Tr \n((p + a-T + im + iF(P) $ ± | 



1±75 



)F(P)(/> + <7-T + im) 



Here, = p^ + eV^, p^ = id^ and Tr stands for the functional trace, i.e. J <i 4 a:tr c trDtr/ 
and F(k) denotes the quark form factor generated by the fermionic zero modes and has an 
explicit form as follows [12]: 



F(kp) = 2z[I Q (z)K 1 {z 



h(z)K (z) - -I^K^z] 



(44) 



where Jo, h, K , and K\ are the modified Bessel functions, z is defined as z = hp/2. When 
k goes to infinity, this zero-mode form factor F(kp) has the following asymptotic behavior: 



F{kp) 
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{kpf 



(45) 



In fact, the nonlocal effective quark-meson interaction can be understood without relying 
on the instanton vacuum [23| . In those cases, the momentum-dependent quark mass can be 
interpreted as a nonlocal regularization in Euclidean space. Hence, various types of the F(k) 
as a regulator with the regularization parameter A ~ 1/p has been used by different authors. 
For example, the dipole-type F{k) is adopted in the study of the pion wave function jj^ . 
while the Gaussian is employed in Ref. |24|. Hence, we utilize in this work the dipole and 
exponential types of the form factors in addition to the instanton-induced zero-mode form 
factor: 



F(k) 



A 2 /(A 2 + k 2 ) (dipole), 
exp(— k 2 /A 2 ) (Gaussian). 



(46) 



The saddle-point approximation of the integrals in the partition function given in Eq. 
leads to the equation for the effective quark mass, which is in general expressed by the 
functional of the external fields V and T. However, since we finally need the term of order 
0(VT), we can consider the saddle-point approximation without external fields: 



*±,fg = S±,/ ff (0) = MfSfg, X±= X 



2V 



U M f 



(47) 



which satisfies the following saddle-point equation 
d 4 k M f F 2 {k){m f + M f F 2 {k)) 



N 
V 



4N r 



(2tt) 4 k 2 + (m f + M f F 2 {k)y 



(48) 



Equation ()48|) contains the momentum-dependent dynamic quark mass Mf(k) = MfF (k). 
Note that the dynamic quark mass Mf is a decreasing function of the current quark mass 
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m f 



251 2(3 ] with meson loops neglected. At small rrif, we may write Mf = Mq + r yrrif with 
M « 360 MeV and 7 « -2 1 . 

Having derived the partition function in Eq. (J4*3*|) . we are now ready for the calculation of 
the magnetic susceptibility Xf horn the instanton vacuum within our framework. In order 
to evaluate x/> we have to compute the following correlation function via the functional 
derivative with respect to the external tensor field: 

5 



\nZ[V,T,m] 



= (0\ip + ^\0) F = Tr(Sa^) (49) 

T=0 



with 



^ P + im + iM(P)' ^ 

We shall consider here only the leading order in the large N c expansion. Using the Schwinger 
method |2lL I22I ] , we solve the following traces: 

Tr 



P + i(m + M(P)) P + irh, 
, / (m + M(p)) rh \ a ■ F 

' (p 2 + (m + M(p)) 2 ) 2 ~ (p 2 + m 2 ) 2 J ~^r GiXV 

+ *{i?T^mw« m( »4 (51) 

where a ■ F = cr^F^. We need to calculate the commutator [p, M(P)], keeping only terms 
to order 0(V). It can be easily computed for the dipole-type form factor. Otherwise, we 
can derive them numerically. Finally, we arrive at the following expression for the magnetic 
susceptibility: 



d A p I mf + Mf(p) nif 



(2vr) 4 I (p 2 + (m f + M f {p)) 2 ) 2 (f + m)) 2 



f d*p M s F{p)F\p)p 
C J (27r) 4 (p 2 + (m / + M / (p)) 2 ) 2 ' 1 1 

where the quark condensate in the chiral limit plays a role of the normalization and is given 
by: 

The infrared region of the first integral in Eq. (j52j) brings out the contribution of order 
Oiviif In rrif) which is almost model-independent. In order to calculate the magnetic sus- 
ceptibility explicitly, we first consider two different regions while integrating: < p < 1 GeV 
and 1 GeV < p < 00. Then we are able to evaluate a part of the free quark analytically. The 



Our preliminary estimate indicates that the contribution of the meson loops to the dynamic quark mass 
is not at all small due to the chiral log corrections which are of order -A- mj In mj . The meson- loop 
corrections are even more important than those of order C(m/). 
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FIG. 1: The nif dependence of the magnetic susceptibility with the form factor of the zero-mode 
type (Eq. 1)44^ ) used. The long-dashed curve draws the contribution of the local vector current, 
while the short-dashed one depicts that of the nonlocal one. The solid curve represents the total 
magnetic susceptibility as a function of the current quark mass. 



second integral in Eq. (j52j) is related to the nonlocal contribution without which the vector 
current is not conserved, i.e. the Ward-Takahashi identity is broken. It arises from the 
nonlocal quark-quark interaction and is called nonlocal current. Summing all contributions, 
we obtain the final result for the magnetic susceptibility of the QCD vacuum at q 2 = 0: 



with 



Xf = Xf + X/° ( 54 ) 

X l r\^Wf)o = (29 + 0.103m f + In ^ ^ [MeV], (55) 

A/ ^/^ /0 ^ * 2vr 2 1000 MeVy L J ' K ' 

X7 nlocal (z^/) = (18 - 0.095771/) [MeV]. (56) 

In principle, we can also derive the magnetic susceptibility at finite q 2 , i.e. the form factor 
of the magnetic susceptibility of the QCD vacuum. 

5. Figure [U shows the magnetic susceptibility x/(^/4V'/)o as a function of the current 
quark mass, with the zero-mode form factor in Eq. (|44|) employed. It is a monotonically 
decreasing function as m/ increases, since the third term with m/lnmj in Eq. ()55|) suppresses 
the first and second ones noticeably as nif gets larger. In fact, its value is around —4.0 MeV 
for up and down quark masses m u ^ ~ 5 MeV, whereas for the strange quark mass m s ~ 200 
MeV it turns out to be —49 MeV. Thus, it almost cancels the leading term in the case of 
the strange magnetic susceptibility Xs- Moreover, the nonlocal part is responsible for about 
40% of the total magnetic susceptibility in the chiral limit and diminishes linearly but slowly 
as mf increases. It implies that it is essential to treat the vector current properly in the 
presence of the nonlocal interaction, in particular, in the case of the magnetic susceptibility. 
As a result, the magnetic susceptibility for the up and down quarks is: Xu,d = 40 ~ 45 MeV, 
whereas for the strange quark we obtain Xs = 6 ~ 10 MeV. As for the up and down quarks, 
the present result is comparable to that of Ref. 0. 0. [J. 
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FIG. 2: A comparison between different types of the form factors. The solid curve is drawn with 
the zero- mode form factor (|44[). while the long-dashed one depicts the dipole-type form factor. The 
short-dashed one is for the Gaussian form factor. 

A comparison between different types of the form factors is depicted in Fig. As already 
discussed in Refs. (13, HiJ , the exponential type of the form factor yields the smallest 
result as usual. It is interesting to see that the dipole type presents a very similar result for 
the magnetic susceptibility to that with the zero- mode form factor (|44)1. It is also found that 
the dependence of the magnetic susceptibility on a type of the form factors is not sensitive. 

6. In the present work, we have investigated the magnetic susceptibility of the QCD 
vacuum within the framework of the nonlocal chiral quark model from the instanton vacuum. 
We first have constructed the effective partition function in the presence of the external 
vector and tensor fields as well as the current quark mass. It was shown that for the constant 
electromagnetic field strength the present result is independent of the path which is selected 
for the gauge connection. We found that the result is a smoothly decreasing function of mf. 
For the up and down quarks m u ^ d ~ 5 MeV we have X-m(^1 d^u,d)o — 40 ~ 45 MeV and for 
the strange quark we obtain Xs^l^sjo = 6 ~ 10 MeV. As for the up and down quarks the 
present results are comparable to the estimate of the vector dominance and QCD sum rule 
approach @, 0, 0| . 
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